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O ! Abstract 



We investigate the routes by which a bath composed of a finite number of oscillators at zero 
temperature approaches the induction of dissipation when it nears the usual limit of dense spectrum 



spread in an infinite interval. It is shown that, when this limit is taken, different distributions of 

& . 

environment frequencies can lead to the same irreversible evolution. However, when we move away 
from it, the dynamics departs from irreversibility in qualitatively different manners. 
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I. INTRODUCTION 



The description of macroscopic irreversible phenomena from microscopic reversible laws, 
in both quantum and classical domains, remains as one of the main challenges of the statis- 
tical physics. In order to take into account dynamical processes involving systems far from 
the equilibrium, like thermalization, transport phenomena, and so on, one must resort to 
kinetic equations m general, kinetic equations govern the effective dynamics of a system 
of interest (a single particle, for example) interacting with a huge number of elementary 
external degrees of freedom of the environment where it is immersed. A kinetic equation 
is obtained from the particular microscopic model describing the interaction between the 
system of interest and its external partners. In classical domain, one usually assumes that 
the dynamics of the global system follows a Hamiltonian equation of motion whereas in 
quantum domain it is governed by the Schrodinger equation. Unfortunately, due to the 
enormous number of degrees of freedom to be considered, usually the exact treatment of the 
microscopic equations of motion is prohibitive and the use of approximations is necessary. 
Thus, the correct understanding of the limits of validity of a given kinetic equation depends 
on the physical conditions imposed to derive it. 

As we remarked above, in the derivation of a quantum kinetic equation it is assumed 
the unitary evolution of the global system composed by the system of interest and the 
environmental degrees of freedom. If a quantum system has discrete eigenenergies, it is 
expected that the evolved state will be found arbitrarily close of the initial state This 
result is the quantum version of the Poincare recurrence theorem. Besides, it is expected 
that the coarse-grained entropy of isolated, finite but large ensembles of quantum interacting 
particles exhibits a periodic or quasi-periodic behavior In general, this last statement 
is valid for quantum systems with integrable classical analogues. So, for a quantum system 
obeying these constraints there is a kind of quantum recurrence time which is so large that 
the larger is such a system. Further, after the seminal work of Anderson [4], it is well- 
known the existence of systems where diffusion is absent by virtue of the inhomogeneous 
distribution of the eigenenergies. 

The last comments illustrate the problem of the extension of the if -theorem, established 
in classical statistical physics by the Boltzmann's pioneer works, to the quantum domain 
(see, e.g., Ref. SB , |6|). The quantum if-theorem intends to set up the conditions to 
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the approach to thermal equilibrium by a macroscopic quantum system and consequently 
irreversible processes involving such a system. In this work, we address the problem of the 
achievement of the irreversible effective dynamics of a system of interest, described by a 
quantum Markovian master equation - an important kind of quantum kinetic equation - 
from the global unitary (and, therefore, reversible) dynamics. Quantum master equations 
are largely applied in several problems in physics, like quantum transport jl-9] and atom- 



photon interactions 



10 



11] . Further, since the seventies, the interest by them has grown by 



virtue of the decoherence program, a theoretical framework devised to give a suitable answer 
to the quantum measurement problem ascribing the disappearance of macroscopic^ quantum 



correlations to the inevitable coupling of a quantum system to its environment 13 ]. 

Since the van Hove's work jl4j . derivations of quantum master equations from first prin- 
ciples have shared a particular ingredient: the number of external degrees of freedom is 
considered infinite, a limit that guarantees the irreversibility of the dynamics of the system 
of interest. In general, one considers the environment modeled by a thermal reservoir pre- 

n 

pared in an initial state that characterizes thermal equilibrium. As discussed in Ref [1J, due 
to the complexity of the knotty action of the environment on the evolution of the system 
of interest, the dynamics of the last can be considered a non- Markovian stochastic process. 
In fact, the evolved state of the global system is completely determined by its initial state 
by virtue of the deterministic nature of the Schrodinger equation. Therefore, the present 
state of the system of interest depends on the whole previous history of its evolution. Under 
suitable limits, it is possible obtain a valid approximation where the memory effects vanish. 
Several methods are employed to obtain a quantum master equation from the particular 



151 ] . Zwanzig 



microscopic model and we cite the projection method developed by Nakagima 

In 



[1 61 ] , and Mori |17| , the cummulant expansion |18fj , and the Feynman- Vernon influence func- 
tional 

The weak coupling limit seems to be ubiquitous in the derivation of a quantum Markovian 
master equation from a microscopic model and this approximation implies to consider the 
global state uncorrelated during all time. However, there is not a uniform way to implement 
the necessary approximations in order to obtain a quantum master equation. In other words, 
the physical conditions that give support to the approximations differ from one particular 
derivation for another. This fact leads some authors to investigate the conditions that a 
proper Markovian quantum master equation could satisfy. These conditions concern, for ex- 
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ample, the preservation of the Robertson-Schrodinger uncertainty relations, the preservation 
of the posit i vity , the existence of a Gibbsian stationary state, and the condition of detailed 
balance 



21 



231 ] . Despite of this, Markovian quantum master equations are successfully ap- 



plied in the description of several phenomena and experiments, especially in quantum optics 
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12). The dynamics of a damped mode of the electromagnetic field, the decay of an 



excited atom by the emission of photons, the process of parametric amplification are few 
examples of processes and phenomena in quantum optics described by Markovian master 
equations. For a dissipative single mode of the electromagnetic field, the system of interest 
is modeled by a harmonic oscillator, in the case of an excited atom decaying to the ground 
state by the emission of one photon, the system of interest is modeled by a two-level system. 
In both cases, the environment is modeled by a bath of harmonic oscillators in thermal 
equilibrium. As a result of the Markovian approximation, the models of damped mode and 
unstable atom predict exponential decay of the excitation of the system of interest, which is 
closely related with the Weisskopf-Wigner theory of radiation decay 24J. This behavior of 
the probability of survival of the excitation as function of time yields a Lorentz-type spectral 
line that was introduced by Breit and Wigner 25[ in order to describe the broadening of 
resonance lines in particle scattering. Despite of fitting properly the experimental data, the 
exponential decay is a byproduct of the Markovian approximation, i.e., the exponential law 
only approximately describes the temporal dependence of the probability of survival of the 
excitation in the system of interest. In fact, one expects that the decay of an excitation in 
quantum mechanics depends quadratically on the time for very short times and is governed 
by a power- law for long times 

The remarks presented above naturally bring some questions to the fore. The van Hove's 
solution to the problem of irreversibility makes use of the limit of infinite number of ex- 
ternal degrees of freedom interacting with a quantum system. However, despite the huge 
number of elementary quantum systems that form a macroscopic object, this number is 
finite. Therefore, if the system of interest is linearly coupled to its external partners, some 
kind of recurrence may happen even though the time necessary for this is possibly greater 
than any relevant characteristic time involved. On the other hand, the region of recurrence 
could be accessible in experiments involving mesoscopic systems. In this work, we intend 
to gain some insight on the achievement of the effective irreversible dynamics of a system 
of interest coupled to external degrees of freedom in the quantum domain. In order to do 
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this, we chose as a laboratory the model of a quantum harmonic oscillator linearly coupled 
to a finite bath of harmonic oscillators in the rotating wave approximation (RWA). These 
external oscillators are prepared in the vacuum state, simulating a bath at null tempera- 
ture, and the main oscillator is prepared in a coherent state. For this case, we can identify 
the conditions on the number of external oscillators, the distribution of their frequencies 
along a finite bandwidth around the frequency of the main oscillator, and the strength of 
the coupling constants in order to obtain a behavior that resembles dissipation within an 
intermediate range of time. All results were obtained numerically in the same sense that 



was done in Ref. [28[ with some important differences. Here, we are interested in studying 
the conditions that lead to an irreversible dynamics; there, the authors study the limits 
of validity of the assumptions generally employed in the derivation of a quantum master 
equation in the Markovian limit. The Markovian limit is also discussed here by the fitting 
of the exponential decay of the energy stored in the main oscillator. 

This work is organized as follows: in the next section, we present the model and determine 
the dynamics of the relevant quantities. In Section III, we study how the relevant parameters, 
like strength of the coupling constants, number of external oscillators and their distribution 
along a finite bandwidth must be adjusted in order to produce an effective dynamics that 
resembles irreversibility. Section IV is reserved to conclusions and final comments. 



II. MODEL AND DYNAMICS 



The model is composed of an oscillator 0\ linearly coupled to a set of oscillators Oj (j = 2 
to N). The oscillator 0\ corresponds to the system of interest and the other ones play the 
role of the environment. The Hamiltonian is 

N N 

H = H + H int , H = y^hv k ala k , H int = (a[a k + d^ij , (1) 



k=l 



k=2 



where a\ and d^ are creation and annihilation bosonic operators, respectively, and v k and 
g k are real coefficients. By displaying this Hamiltonian in the matrix form 
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H = h I a\ a 2 



A' 



vi g 2 ■■■ g N 

g 2 v 2 

: '•• ; 

9n ■ ■ • v N 



«2 



\ 0-N J 



(2) 



and denning 



/ 



M 



Ml,! .U; , ••■ M 1>N 

M 2 i M 2 o ■■■ M 2 n 



\ 



(3) 



\ M^i .1/y, • • • M N>N J 
real orthogonal matrix satisfying 



f vx g 2 ■■■ g N ^ 



M 



/ 



M 



(4) 



g 2 v 2 

: '•. ; 

y 9n • • • v N J y 

it becomes easy to see that the normal modes operators are obtained through the variables 
transformation 



wi ••■ 

u 2 

: '•• : 

■■■ u N 



\ 



J 



A 2 



\ An J 



M 1 



fa ^ 



a 2 



(5) 



since Hamiltonian (CQ) can be written as 

/ 



H = h(A\ A\ ■■■ A 



ui ••• 

u 2 

: '•• : 

••• u N 



\ 



fA\ 

A 9 



\ An J 



N 



(6) 



k=l 



The orthogonality of M can be used to show that the operators A k and A^ obey the usual 
commutation relations for bosons. The vacuum of the normal modes is the same as that of 
the original modes, what is easily proved by noting that the annihilation operators of the 
normal modes are linear combinations of the annihilation operators of the original modes. 
Let us assume that the initial state of the whole system is given by 



|V (0)} 



«i, a 2 , 



(7) 



where the vector state |oti, a 2 , q;jv) indicates that oscillator Ok is in the coherent state 
with amplitude ak- Using Eq. (JHJ), we see that the coherent states related to the original 
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operators a k are also coherent states associated to the normal modes operators A k . 
\a x ,a 2 , ...,a N ) = exp (a k a{ - a* k a k ^j |0, 0, 0) 



A? 



= exp r£ (A k Al-A* k AM |0, 0> 

= |Ai, A 2 , Ajv) , (8) 
where |Ai, A 2 , Ajv) is a coherent state related to the normal modes with amplitudes 

N 

A k = J2 a oM jk - (9) 

Thus, 

\i/)(t)) = \A 1 e- lu ' lt ,A 2 e~ iul2t ,...,A N e- iulNt ) 

= exp HT (A fc e-^4-A*e^)j |0, 0, 0> 



exp (^>2 a k (*) 4 _ a l (*) fl fe^ 1°) °) 0) 
|«i (*) ,ot 2 it) , ...,ajv (*)) , 



(10) 



N N 



-iujit 



(11) 



where 

a k (t) = J2J2 a i (°) M i>* M M e ~ 
j'=i /=i 

Since all oscillators remain in coherent states, the environment does not produce decoherence 
on the system. If the central oscillator starts in the coherent state with a% = 1 and the 
environment is in vacuum state, the energy of oscillator 0\ evolves as 

JV 



E\ (t) = (%v\a\ai) (t) = hvi \azi (t)| z = %v\ 



(12) 



As can be shown by using results in Ref. |29|, expression (fT2l is also valid for the central 
oscillator starting in the Fock state with one excitation and the environment in vacuum. 
Unlike the coherent state, for such initial condition the system loses purity. 



III. THE WAY TO IRREVERSIBILITY 



From now on, the system of units is chosen so as %V\ = 1, and the values of all physical 
quantities below are given with respect to this system. If the environment is composed of 
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a set of oscillators with equally spaced frequencies varying in the range from 1 + p/2 to 
1 — p/2, and the coupling constants are equal, decaying with 1/yN — 1 when the number 
of environmental modes grows, we get, for even N, 
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(13) 



where /i = p/ (N — 2). The form of is analogous to the one assumed in master equation 
derivations in order to achieve finite environmental effects when the thermodynamic limit 
N — > oo is taken (see, for example, Ref. 



The dynamics of Ei (t) depends on the distributions of M\ k and Uk, which, although 
usually calculated only numerically, can be understood by considering the general situation 
as intermediate between analytically tractable cases. In Figs. ([Q), we exemplify the change- 
over of the set |M^ fc }. For small g/p, the matrix in expression ( fl3l) approximates a diagonal 



matrix with a two fold degenerate eigenvalue 1, leading to M\ 



N/2 



M 2 

1V1 l,N/2+l 



1/2 and 



M 2 k = otherwise. For large g/p, the set |M^ fc } becomes like the one for p = 0, where the 
only nonzero elements are M\ x = Mf N = 1/2; this can be understood by noticing that if M 
diagonalizes H then it also diagonalizes H/g, whose diagonal elements approach each other 
when g/p increases. When the ratio g/p grows from zero, {Mf k } varies from the former 
situation to the latest. 

Let us start by investigating the intermediate case g/p = 0.2. For t = 0, we have E\ = 1. 
From the mathematical point of view, this result is a consequence of the orthonormality of 
M. As time progresses, the terms in the summation in expression (fl2l) gain different phases, 
leading to the decreasing of E\ it). As shown in the Appendix, the time scale for this decay 
is given by 

r d = ± - 3.232, (14) 
8 



where 7 and (3 are parameters obtained through Lorentzian and linear fittings on the points 
of Figs. (E^,) and ([2b), respectively. Notice that this scale is in agreement with the evolution 
in Fig. ()2t). A partial revival occurs around 

r, - y (15) 

(see Fig. ([2H)). In order to understand this revival, we start by observing that the distri- 
bution in Fig. (|2b) is almost linear. Accordingly, the differences Uj — Wj-i for j = 2 to 
N are, in good approximation, given by their mean value (con — u)i) / (N — 1), and, when 
t = T r , e~ luJjt and e -1 ^ -1 ' accumulates a phase difference close to 2tt; then every terms in 
the right hand side of Eq. ( fl2l) will be approximately in phase. Other revivals occur for 
t = nr r , where n is an integer greater than 1; when n increases, they go spreading. Of course, 
revivals are not expected for a model of a dissipative environment. In fact, they become 
unimportant when the density of environmental modes tends to infinite, since in this limit 
r r — > 00. Another way to avoid these relatively localized revivals is to randomly distribute 
the frequencies of the environmental modes, that leads to permanent fluctuations that may 
be associated to a spreading of the revival. As may be seen in Figs. ([3^) and ([3b), these 
fluctuations decrease as N increases, becoming unimportant in the thermodynamic limit. 

In Fig. (jUi), we display the evolution of Ei (t) for g/p = 0.8, which may be regarded a 
large value for this ratio. The oscillations of the energy between system and environment 
are substantial. According to Figs. (Hti), (Hb) and (TJb), these are well approximated by 
preserving only the most relevant terms in Eq. (fl2l) : 

E 1 [t) w jM^e-^ 1 ' + Ml N e~ iulNt \ 2 w 0.367 [1 + cos (1.714*)] . (16) 

Such oscillations, which are not expected for a model of dissipative environment, are re- 
lated to the finite range of environmental frequencies, decreasing when p is enhanced, as is 
exemplified by comparing Figs. (jUi) and (JSJ). This finitude becomes important when the 
coupling strength with the environment increases. 

The case of weak coupling is treated in Figs. ([6b.) and ([6b), prepared with g/p = 0.01. 
Now, 

Ex (t) w \Ml N/2 e~ iujN ^ + Ml 1+N/2 e- luJl+N / 2t \ 2 w 0.469 (1 + cos [(0.002*)]) , (17) 

agreeing with Fig. (Et). The inflection in Fig. ()6b) indicates that the difference between 00^/2 
and oj\ + n/2 is shorter than the other ujj — ujj^i differences, what reflects the approximation 
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to the condition g = 0, where there is a degenerate eigenvalue: ujn/2 = = 1- The 

freezing of the dynamics associated with the decrease of g comes from the increasing of the 
period of the main oscillation. It is noteworthy that, if the resonant environmental mode is 
taken out, the dynamics for low g/p will be qualitatively different (see Figs. (J7b) and (EH)). 
As g/p goes to zero, the eigenvalues con/ 2 and Ui + n/ 2 remain non-degenerate and M^r N+1 y 2 
approaches the unit. With the help of Figs. (JTK) and ([7b), we see that 

-^1 \t) ~ | iW l,(7V+l)/2-l e + iW l,(Af+l)/2 e + iW l,(W+l)/2+l e k " | 

« (0.968 + 0.019 cos (O.OlOt)) 2 , (18) 

where we used M 1 2 (7V+1)/2 _ 1 « M 2 (Ar+1)/2+1 and£«;(iv + i)/2-W(jv+i)/2-i ~ ^(iv+i)/2+i-^(iv+i)/2- 
The standstill of i?i (t) is achieved through the increasing of Mf, N+1 y 2 . 



IV. CONCLUSIONS 



Irreversibility is often expected for the dynamics of a system affected by the environ- 
ment. This is usually achieved by considering the set of environmental degrees of freedom 
in certain limits, as dense spectrum spread in an infinite interval. In order to investigate 
the way the system's dynamics approaches the irreversibility, the environments we analyzed 
here are finite. As we show, the expected behaviors may be reached through different ways. 
The dynamics we calculated depends crucially on the sets of eigenvalues and eigenvectors 
of the Hamiltonian of system plus environment. Although these sets are not obtainable 
analytically, we analyzed their structures by considering them as varying between two ana- 
lytically tractable cases. This way, the numerically calculated dynamics can be more deeply 
understood. 

Since the bath discussed is at null temperature, we expect that it produces dissipation. 
In order to obtain such a dynamics, the relation between the interval where the environment 
frequencies are distributed and the strength of the coupling to the environment must be 
properly chosen. If the coupling dominates, the environment approaches the situation where 
all oscillators have the same frequency, leading to a collective behavior of the environment 
oscillators that induces an evolution similar to the one concerning a single oscillator highly 
coupled to the system, i.e., fast periodic energy flow between system and environment. In an 
suited situation, the relation between coupling strength and bandwidth must be such that the 
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effects of the oscillators out of this bandwidth are negligible, due to detuning. Approaching 
this case, the fast oscillations described above progressively decrease, until we get a dynamics 
that resembles dissipation. The role of a large (in the limit, infinite) bandwidth in the model 
is to avoid the collective behavior of the oscillators of the environment. 

Since the bath is finite, the energy that flows from the system returns, in part, in finite 
times. The form of this flow back depends on the distribution of the bath modes: for equally 
spaced frequencies, there are well defined partial revivals; for randomly distributed frequen- 
cies, there are no revivals, but the residual oscillations are relatively large. As we approach 
the condition of dense spectrum, irreversibility is achieved, for equally spaced frequencies, 
through the growing of the revival time, while, for randomly distributed frequencies, through 
the decrease of the residual oscillations. 

For relatively strong coupling, the importance of each single bath oscillator for the system 
dynamics is small. Nevertheless, for weak coupling, the presence of one single resonant 
mode is crucial for the definition of the shape of the dynamics. The portion of energy 
that two coupled linear oscillators can exchange depends on the relation between coupling 
and detuning. For weak coupling, this portion will be small for non-resonant oscillators, 
while the whole energy is aways exchangeable for the resonant one. Accordingly, if there is a 
resonant mode, the interaction of the system with this mode dominates. Thus, the dynamics 
approaches the one where the system interacts with a single loosely coupled oscillator, leading 
to a large period for the flow of energy. This period increases with the number of oscillators, 
paving the way for irreversibility. Already for absent resonant mode, only a small part of 
the energy can flow between the system and the environment. Exact resonance is hardly 
expected in a real case, but if there is one oscillator for which the ratio between dephasing 
and coupling is low and much lower than for the other, the dynamics may be close to the 
one with the resonant mode. 

The limit conditions leading to irreversibility can be reached by modeling the bath with 
different frequency distributions. Although the dynamics induced are the same when the 
limit conditions are assumed, if we move away from them the evolution of the system can 
vary qualitatively. 
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Appendix A 



According to Eq. (112]) . the energy of oscillator 0\ may be written as 



E 1 (t)=hv 1 \a(t)\ 2 , (Al) 

where 

a (t) = J^M^e-****. (A2) 

3=1 

In this Appendix, we employ the concept of Fourier transform of a function / (t), 

+00 

F(k) = -J=Jf (t) e M dt, (A3) 

—00 

and its inverse, 

+00 

1 



f(t) = -= I F (k) e-^dk, (A4) 



-00 



to analyze o it). We will explore possible relations between eqs. 0A2p and (IA4[) . 

For the parameters of Fig. (j2J), it may be assumed that Uj ~ u + f3j, where ujq ~ 0.4902 
and (3 ~ 9.899 x 10 -3 are obtained through linear regression. Using this approximation, we 
can write 

a (t) « e~ iu:oi J2 M l (j) e ~ mt i ( A5 ) 
i=i 

where Mi (x) is a continuous function defined for every real number x such that M\ (j) = 
Mi j for j = 1, 2, • • • , N. By changing the summation into an integral and considering that 

Mf (x) is negligible for x < or x > N, we get 

+00 

a (t) « e~ iaJot / M 2 (z) e" it/3a; dx, (A6) 



which, after the variables transformation k = (3x, leads to 

v'^ , ,2 fk\ ,,,, 



e i0J °V (t) « -= / -^M 2 - ) e^dife. (A7) 
\/27r j p VP, 
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$ (k) = ^-M? ( \ ) (A8) 



By comparing eqs. (1A4I) and (1A7I) . we see that 

7" *" 1 V/3 

is nearly the Fourier transform of 

(t) = e iwot a (t) . (A9) 

Of course, we can substitute a (t) for (£) in Eq. ( 1A1I) without any modification in i?i (£). 
Now we assume that Ml (x) is approximately the Lorentzian function 

Ml (x) « ^— (A10) 

where A m 0.3886, 7 m 15.63 and x ~ 51.50 are obtained by fitting the data in Fig. (J2J). 
Thus, 

(A: - f3x Q ) 2 + {hf 

what is the Fourier transform of 



*(*)« „ :,r'^ 2 ; (aid 



w W^e"^ 1 * 1 , (A12) 

leading to 

Ei(t) ^ /ivrV^e- 27 ^ 1 ' 1 . (A13) 
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(c) (d) 

FIG. 2: (a) and (b) give the sets fc j and {uj^}, respectively, (c) and (d) show, in different 
time scales, the evolution of E\. For all figures, H is given by Eq. (fT3|) with N = 102, p = 1 and 
g = 0.2. 
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FIG. 3: Ei as a function of time for hv\ = 1, TiVk randomly distributed in the interval [0.5,1.5], 
g k = 0.2/y/N - 1 (k = 2 to N) and (a) N = 102, (b) N = 1002. 
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(c) (d) 

FIG. 4: (a) and (b) give the sets |Af^ fc | and {uik}, respectively, (c) and (d) show, in different 
time scales, the evolution of E\. For all figures, H is given by Eq. (|13p with N = 102, p = 1 and 
g = 0.8. 
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FIG. 5: E\ as a function of time for H given by Eq. (|13p with N = 102, p = 2 and g = 0.8. 
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FIG. 6: (a) and (b) give the sets fc j and {wfc}, respectively, (c) shows the evolution of E\. 
For all figures, H is given by Eq. (fT3|) with N = 102, p = 1 and g = 0.01. 
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FIG. 7: (a) and (b) give the sets j^f^ fc j and {uJk], respectively, (c) and (d) show, in different 
time scales, the evolution of E\. For all figures, H has the form given in Eq. f 1 1 3 [) . except for 
second line and second column absent. Here, N = 101, p = 1 and g = 0.01. 
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